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AKSTIUCT 


A  hypothesis  concerning  the  ratio  of  the  exchange  coefficients 
tor  heat  and  momentum  it  formulated.  From  the  basic  assumptions  of 

W 

the  Mon  in-Obukhov  similarity  theory  and  the  theory  of  free  convec- 
tioM,  it  is  shown  that  the  ratio  of  the  coefficieBts  is  a  two-part 
function  of  the  nondimens ional  logarithmic  wind  shear  and  the  gradi¬ 
ent  Richardson  number.  Experimental  data  tend  to  corroborate  the 
theoretical  values  derived  in  this  study  within  the  limits  of 
experimental  error. 
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INTROMCTKA 


Contemporary  hypotheses  pertaining  to  the  turbulent  processes  of  the 
surface  boundary  layer  are  based  upon  the  principles  of  dynamic  similarity  of 
flows.  True  dynamic  similarity  requires  that  the  measured  quantities  of  one 
system  go  identically  or  by  an  exact  ratio  into  those  of  another,  implying  a 
constant  ratio  of  forces  in  the  two  systems.  For  the  surface  boundary  layer, 
this  requires  that  the  wind  and  temperature  profiles  be  similar,  a*  well  as 
the  vertical  fluxes  of  heat,  momentum,  and  water  vapor. 

The  structure  of  turbulence  in  the  boundary  layer  is  dependent  upon  the 
thermal  statification.  Richardson  (1920,  1925)  found  that  the  generation  or 
suppression  of  turh.iloni**  could  be  expressed  as  the  ratio  of  work  done  by  buoy¬ 
ancy  forces  against  gravity  forces  where  turbulence  increased  as  siaLliity  de¬ 
creased.  The  converse  was  also  found  to  be  true.  Richardson’s  criterion,  as 
originally  derived,  assumed  that  the  exchange  coefficients  for  heat,  momentum, 
and  water  vapor  were  equal,  allowing  the  Richardson  number  to  be  expressed  in 
generalized  form.  Current  thinking  indicates  that  the  exchange  coefficients 
are  not  equal,  so  that  the  Richardson  number  may  be  represented  by  either  its 
flux  or  gradient  form.  The  two  Richardson  numbers  are  related  by 

K„ 

Rf  =  Ri  -li  (1) 

KM 

where  Rf  is  tiie  flux  form,  Ri  is  the  gradient  form  and  K||/Kn  Is  the  ratio  of 
the  exchange  coefficients  for  heat  and  momentum. 

Batchelor  (1953)  demonstrated  th8t  the  Richardsoa  number  possessed  the 
characteristics  of  true  dynamic  similarity,  which  allows  it  to  be  used  as  a 
basic  parameter  for  the  investigation  of  turbulent  processes  in  the  surface 
boundary  layer.  This  leads  to  the  premise  that  the  ratio  of  the  exchange  co¬ 
efficients  is  also  a  similarity  ratio  of  major  importance  owing  to  the  nature 
of  its  relationship  to  the  Richardson  numbers.  If  dynamic  similarity  of  flows 
indeed  exists  in  the  planetary  boundary  layer,  it  should  be  possible  from  the 
basic  definitions  to  write  an  expression  for  the  ratio  of  the  exchange  coef¬ 
ficients  for  heat  and  momentum. 

The  purposes  of  this  paper  then,  are  to:  (1)  examine  the  role  that  K^/K^ 
assumes  in  the  analysis  of  wind  and  temperature  profiles  observed  in  the  planet¬ 
ary  boundary  layer,  (2)  examine  the  dependency  of  diabatie  boundary  layer  hy¬ 
potheses  upon  Kq/K^,  and  (3)  to  present  an  expression  for  computing  the  ratio 
in  terms  of  similarity  parameters.  In  addition,  values  of  Kq/KM  as  determined 
from  experimental  data  are  presented. 

SIMILARITY  CONCEPTS 

The  application  of  dynamic  similarity  concepts  ta  surface  boundary  layer 
flow  is  usually  attributed  to  Monin  and  Obukhov  (1954),  even  though  the  earlier 
work  of  Sverdrup,  Rossby  and  Montgomery,  and  Holtzman  (see  Deacon,  1955  and 
Lettau,  1949)  was  based  upon  this  postulate  of  aerodynamics.  The  conjugate 
law  for  the  wind  profile  in  the  surface  boundary  layer  may  be  written  as 
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*1  =  «a  (s>  <2) 

*7  k7 

where  V  is  the  wind  speed,  /  is  height,  u*  is  the  friction  velocity,  k  is  liar- 
man's  constant  and  S  is  a  stability  correction  factor  known  as  the  diahatic 
influence  function.  Monin  and  Obukhov  (195-1)  determined  that  dynamic  siMilar- 
i tv  of  flows  was  controlled  by  the  unique  existence  near  the  ground  of  the 
friction  velocity  »*,  a  scaling  temperature  T*  and  a  scaling  length  L.  All  the 
equivalent  measured  values  could  then  be  expressed  as  functions  of  these  param¬ 
eters.  With  the  aid  of  the  dimensional  constants  g/Q  and  H/cpp,  L  and  T*  were 
defined  as 


and 


L 


ul  V9 

kgll 


•ye  -  ,  1  H 

k\  CPP 


(3) 


(4) 


where  g  is  the  acceleration  of  gravity,  0  is  potential  temperature,  H  the  ver¬ 
tical  heat  flux,  Cp  the  specific  heat  of  air  at  constant  pressure,  and  p  the 
mean  density. 


From  the  above  definitions,  an  exact  statement  describing  the  shape  of 
the  wind  and  temperature  profiles  in  a  diabatic  surface  boundary  layer  ray  b«_ 
written  as 


(3) 


where  z/L  is  a  scaling  ratio,  and  the  flux  Richardson  number  is  given  by 

■  t 

i 


Cp0r  (dv/dz) 


(6) 


The  wind  profile  may  then  be  determined  by  integration  to  be 


where  z0  is  a  constant  of  integration  known  as  the  roughness  length  and  a  is  a 
universal  constant. 
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Since  the  expressions  for  both  L  and  Rf  contain  the  heat  flux  H,  a  param¬ 
eter  difficult  to  obtain  by  inference  or  direct  measurement,  Eq.  (5)  may  be 
restated 


I*  -  S  Hi 


(8) 


where  Ri  is  the  gradient  Richardson  number  defined  as 


(9) 


and  L'  is  e  gradient  length  (Panofsky,  Blackadar,  and  McVehil,  i960)  given  by 


L' 


dv 

Vk  8  5? 


z 


d6 


(10) 


and  related  to  L  by  L* 


-  L  kH/KM. 


The  wind  profile  may  now  be  stated  as 


(11) 


where^  is  a  universal  function.  A  similar  expression  may  be  written  for  the 
temperature  profile. 

Thus,  from  easily  determined  parameters  which  are  functions  of  the  verti¬ 
cal  gradients  of  wind  and  temperature,  the  turbulent  processes  in  the  surface 
boundary  layer  may  be  investigated.  Accurate  evaluation  of  the  vertical  fluxes 
from  profile  measurements  is  entirely  dependent  upon  the  ratio  of  the  exchange 
coefficients . 

Another  factor  influencing  the  dynamic  similarity  of  flows  in  a  diabatic 
surface  boundary  layer,  particularly  as  the  atmosphere  becomes  less  stable,  is 
the  effect  of  convection  upon  the  turbulent  processes.  It  is  generally  accept¬ 
ed  that  there  exists  a  transition  zone  between  the  forced  and  free  convection 
regimes  at  -0.02>  Ri>  -0.05.  In  fully  forced  convection  (Ri  >  -0.02),  Priestley 
(i960)  has  shown  that  the  potential  temperature  gradient  can  be  stated  as 
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(12) 


de  ht  SH 

3*  *  "  cppkM»  2 

Beyond  the  transition,  the  profile  no  longer  follows  the  z“’  law  of  forced  con¬ 
vection,  but,  as  shown  by  Priestley  (1954),  becomes  proportional  to  height  to 
the  minus  four-thirds.  The  potential  temperature  gradient  in  free  convection 
then  becomes 


de 


-X/3  -4/3 

(t:/e)  z 


(13) 


where  h  is  a  constant  whose  value  remains  to  be  determined.  In  the  eva|uation 
of  h,  Pri“«t*«»y  0955)  introduced  a  reduced  or  dimensionless  heat  flux  H  in  the 
form 


allowing  the  flux-gradient  reciprocal  dependence  between  free  and  forced  convec¬ 
tion  to  be  given  by 


»> 

H  =  h 


(free)  (15) 


H  =  k2  |Ri|  l forced)  (16) 

where  II  is  a  constant  with  height  in  the  free  convection  regime  and  a  variable 
in  the  forced  regime. 

Egs.  (11)  through  (16)  were  evaluated  using  Swinbank's  (1955)  date  in 
terms  of  ft  and  Ri.  An  extension  of  this  inquiry  by  Taylor  (1956)  resulted  in 
establishing  a  value  for  h  such  that 

* 

H  -  h  =  0.79  ±  0.04 

According  to  Priestley  (i960),  measurements^of  the  heat  flux  were  some  10  per 
cent  low  and  a  tentative  adjusted  value  of  H~h  =  0.9  was  adopted  for  the  free 
convection  regime. 

Heat  flux  measurements  by  Priestley  (1955)  and  Taylor  (1956)  and  tempera¬ 
ture  profiles  by  Webb  (1958)  indicate  the  height  to  the  minus  four-thirds  law 
for  free  convection  is  valid  for  z/l.  <  0.03  through  a  range  up  to  about  30 
times  this  value.  Egs.  (12)  and  (13)  then  represent  the  behavior  of  the  temper¬ 
ature  structure  in  the  surface  boundary  layer  with  the  scale  height  of  the 
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transition  from  forced  to  free  convection  proportional  to  the  Mon’.i-libukhov 
scaling  length  L. 

The  transition  between  forced  and  free  convection  can  be  considered  smooth 
since  moving  turbulent  elements  can  traverse  the  junction.  With  no  smoothing, 
it  can  be  shown  from  the  diabatlc  similarity  theory  that  a  Junction  height  is 
given  by 


\  -  -*4  h‘a  •  (IT) 

Eq.  (17)  yields  an  arbitrary  junction  height  of  z/L  =  -0.0317  if  it  is  assumed 
that  k  =  0.4,  h  =  0.9  and  Kjj/Kj.  is  set  equal  to  1.  The  value  of  the  gradient 
Richardson  number  is  also  found  to  be  -0.0317  at  the  unsmoothed  Junction. 

THE  RATIO  OF  THE  EXCHANGE 
COEFFICIENTS  FOR  HEAT  AND  MOMENTUM 

Now,  from  the  related  diahatic  similarity  and  free  convection  theories, 
we  can  formulate  an  expression  for  the  exchange  coefficients  for  heat  and 
momentum.  The  exchange  coefficient  hypothesis  states  that  the  mean  flux  per 
unit  area  of  a  conservative  quantity  such  as  heat,  momentum,  or  water  vapor  is 
proportional  to. the  gradient  of  the  mean  value  of  the  quantity,  that  is 


mean  flux  per  unit  area  =  -K 

where  X  is  the  exchange  coefficient,  dE/dn  is  the  gradient  of  the  mean,  and  n 
is  the  direction  normal  to  the  surface.  For  turbulent  flow  in  the  boundary 
layer,  K  is  dependent  upon  time  and  location.  A  differentiation  between  the 
various  coefficients  is  necessary  to  adequately  describe  turbulent  processes  in 
the  boundary  layer.  These  c*n  be  stated  as 
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(20) 


where  K^,  Kfj,  and  are  the  coefficients  for  water  vapor,  heat,  and  momentum, 
respectively. 

Considering  only  heat  and  momentum,  the  exchange  coefficients  can  be  shown 

to  be 
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r  - 

h  30 

3s 

(21) 

and 

r  -  22 

3V 

3* 

(22) 

which  can  be  verified  by  direct  observation.  Since  independent  observations  of 
the  terms  necessary  to  provide  a  solution  for  Eqs .  (19)  through  (22)  are  diffi¬ 
cult,  if  not  sometimes  impossible,  to  obtain  by  direct  measurement,  an  indirect 
determination  of  the  ratio  of  to  Ky  can  be  derived  from  the  basic  framework 

of  the  similarity  and  free  convection  theories. 

From  Eqs.  (1)  and  (9)  it  is  seen  that 

»  .  R1 £  S  2 

f  "m  KM0/| vj 

(23) 

Substitution  of  Eq.  (19)  and  (20)  yields 

R  _ _ H&_ 

c  T0  3V  • 

P  3^ 

(24) 

By  definition  u2  -  x/p,  so  that  Eq.  (19)  may  be  restated 

as 

Rf  “  pc  u.^6  3V 

P 

(25) 

Introducing  Priestley’s  (1955)  reduced  heat  flux  H,  it  is  seen  that 


_  Hrz2(k/9)**  jilZlzl  3/? 


3z 


(26) 
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Multiplying  and  dividing  by  k?  (4V/I*)*,  introducing  the  diabetic  influence 
function  S,  and  rearranging  terms,  Eq.  (26)  becomes 


Rf  *  - 


2  *  W. 

Lll 


k  e  ( jv/ji) 


4§Z>«1 


v2 


(27) 


Simplifying 


Since  by  definition  R^, 


v7 


Rl. 


=  Ri  (Ku/K),  then 

H  M 


(28) 


If  the  junction  between  forced  and  free  convection  is  assumed  to  occur  at  Mi  = 
-0.0317,  then 

(a)  for  forced  convection  where  |Ri]  <0.0317  and  H  =  |Rf  ~^k3 


h 

•Si 


kL  s* 
IriI^ 


s 


V 


(b)  for  free  convection  where  iRil  >0.0317  and  H  =  h  =  0.9 


(29) 


•Si 


0.9  S2 

“kT" 


(30) 


It  is  seen  that  the  ratio  Kj)/K^  may  be  calculated  from  two  of  the  basic  param¬ 
eters  of  similarity  theory,  the  Richardson  number  and  the  diabatic  influence 
function  S. 


NUMERICAL  EVALUATION  OF  THE 
RATIO  OF  THE  EXCHANGE  COEFFICIENTS 

The  ratio  of  the  exchange  coefficients  may  be  evaluated  numerically  by  use 
of  any  of  the  contemporary  hypotheses  for  the  diabatic  boundary  l8yer.  Since 
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none  of  the  dozen  or  so  hypotheses  are  universally  accepted,  the  choice  is  rather 
arbitrary.  Application  of  Eqs.  (29)  and  (30)  is  dependent  upon  the  definition 
of  S  assuaed  as  a  stability  parameter  for  solution  of  the  log-linear  profile. 

By  definition 

,  .  15  W 

dz  (31) 

the  solution  to  which  is  available  only  with  independent  measurements  of  stress. 
The  Mon  in-Obukhov  (1954)  hypothesis  assumed  that 

S  -  (1  +  a*)  (32) 

which  is  a  first  approximation  for  small  z/L  only  and  consequently  is  not  valid 
in  the  free  convection  regime.  Interpolation  formulae  to  span  the  transition 
zone  between  forced  and  free  convection  and  fit  at  both  small  and  large  z/L  have 
been  formulated,  those  of  Businger  (1959,  1961),  Priestley  (1960),  Webb  (1960) 
and  the  KEYPS  function  of  Panofsky  ( 1963)  being  the  most  notable. 

Depending  upon  the  choice  of  an  interpolation  scheme,  slightly  different 
values  of  S  are  obtained  for  the  same  value  of  Ri  or  z/L',  causing  the  stress 
determination  as  obtained  from  Eqs.  (7)  or  (9)  to  differ.  The  same  holds  true 
for  the  determination  of  KH/KH  from  Eqs.  (29)  and  (30).  The  differences  among 
the  expressions  for  the  diabatic  influence  functions  as  proposed  by  various  in¬ 
vestigators  are  currently  unresolved. 

In  this  particular  study,  the  KEYPS  function  was  chosen  to  evaluate  Ku/Km. 
Panofsky  (1963)  has  shown  that  for  a  diabatic  surface  boundary  layer,  the  diaba¬ 
tic  influence  function  from  the  KEYPS  expression  is  given  by 

S*  -y  l  S3  -  1  (33) 

where  y  is  a  universal  function.  In  terms  of  the  gradient  length  L*.  this  may 
be  expressed  as 

S*  -  y'  J,  S3  -  1  (34) 

where  y'  is  another  universal  function.  From  similarity  theory  it  can  be  shown 
that 

z  _ Ri _  . 

L'  "  (1  -  Y  '  Ri)*t 

and 


S  -  (1  -  Y  *  RI)  (36) 
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allowing  determination  of  the  various  profile  parameters  if  the  vertical  gradi¬ 
ents  of  wind  speed  and  temperature  are  measured. 

Numerical  evaluation  of  Kh/Km  was  accomplished  using  high-speed  digital 
computer  techniques.  Eqs.  (29),  (30),  (35),  and  (36)  were  programmed  to  pro¬ 
vide  a  solution  in  terms  of  S,  -Ri,  -z/L'  and  the  equivalent  values  of  Kh/Km  in 
in  the  stability  range  0  >  Ri  >  -10,0.  Kh/Km  as  a  function  of  -Ri  and  -z/L*  is 
shown  in  Figures  1  and  2. 

DISCUSSION 

The  solution  of  many  problems  of  boundary  layer  flow  depends  upon  the  ex¬ 
act  numerical  value  of  the  ratio  of  the  exchange  coefficients  for  heat  and  mo¬ 
mentum.  Early  experimenters  such  as  Richardson  (1920,  1925)  assumed  that  KH/Kg 
=  1,  an  assumption  not  Justified  by  current  hypotheses.  The  evolving  similarity 
theory  of  Monin-Obukhov  utilized  this  assumption.  Extension  of  the  Monin-Obuk- 
hov  hypothesis  by  the  interpolation  formulae  tended  to  indicate  that  Kh/Km  was 
closer  to  1.3  in  the  free  convection  regime. 

Experimental  values  of  the  ratio  of  the  exchange  coefficients  have  been 
published  by  a  number  of  investigators,  notably  Rider  (1954),  Swinbank  (1955), 
Deacon  (1956),  and  Senderikhina  (1961).  Published  values  range  from  1.06  to 
1.67,  although  Lumley  and  Panofsky  (1964)  quote  unpublished  data  by  Priestley 
showing  that  values  as  large  as  3.0  have  been  observed  in  the  boundary  layer. 
Brooks  (1963)  found  values  as  high  as  2.5  during  extremely  unstable  periods  over 
an  irrigated  field,  but  according  to  Dyer  and  Truitt  (1962)  and  Dyer  (1963),  the 
process  in  operation  when  Brook's  data  were  collected  was  heterogeneous,  probably 
accounting  for  the  large  fy/fy  ratio. 

Lettau  (1957)  analyzed  two-dimensional  wind  and  temperature  data  obtained 
during  the  Great  Plains  Turbulence  Field  Program  and  determined  the  ratio  of  the 
exchange  coefficients  from  the  Bowen  ratio  and  a  similarity  assumption  referred 
to  as  the  "Ground-flux  Ratio."  It  was  assumed  that  conditions  valid  for  the 
neutral  case  held  true  for  the  forced  and  free  convection  regimes  also.  This 
led  to  biasing  the  data  to  the  extent  that  computed  values  of  K^/Kg  in  the  range 
Ri  <-0.0317  are  abnormally  high  as  shown  in  Figure  3. 

All  evidence,  both  empirical  and  experimental,  indicate  that  the  ratio  of 
the  exchange  coefficients  for  heat  and  momentum  varies  between  unity  and  1.3 
from  the  neutral  case  to  the  upper  limit  of  the  free  convection  regime  which 
lies  in  the  stability  range  -0.1>Ri>  -1.0  according  to  Nebb  (1956)  and 
Townsend  (1962).  This  is  considered  to  be  a  second  transition  zone  and  separates 
the  free  and  "windless"  or  natural  convection  regimes.  The  temperature  profile 
no  longer  follows  the  minus  one-third  law  of  free  convection  but  obeys  a  z“a 
law  according  to  Lumley  and  Panofsky  (1964).  This  is  borne  out  by  the  failure 
of  any  of  the  diabetic  models  to  predict  stress  and  roughness  lengths  accurately 
at  large  negative  values  of  Ri,  since  the  interpolation  formulae  are  tailored  to 
free  convection  theory  and  the  minus  one-third  potential  temperature  profile  of 
that  stability  regime.  The  effect  on  predicted  or  measured  values  of  KH/Kg  is 
not  known. 
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FIGURE  I.  THE  RELATION  BETWEEN  THE  RATIO  OF  THE  EXCHANGE  COEFFICIENTS  FOR  HEAT 
AND  MOMENTUM  AND  THE  GRADIENT  RICHARDSON  NUMBER. 
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FIGURE  2.  THE  RELATION  BETWEEN  THE  RATIO  OF  THE  EXCHANGE  COEFFICIENTS  FOR  HEAT 
ANO  MOMENTUM  ANO  THE  KEYPS  RATIO. 


FIGURE  3.  RATIO  OF  THE  EXCHANGE  COEFFICIENTS  vt  Rl  AT  1.6  METERS 


CONCLUSIONS 


The  empirical  solution  for  the  ratio  of  the  exchange  coefficients  for  heat 
and  aomentum  yields  numerical  values  in  good  agreement  with  the  published  results 
derived  from  experimental  data.  It  must  be  emphasized  that  the  solution  is  val¬ 
id  only  for  a  homogeneous  and  stationary  diabetic  boundary  layer.  Other  assump¬ 
tions  used  for  the  solution  are  that:  (1)  the  basic  similarity  and  free  convec¬ 
tion  theories  are  valid;  (2)  that  the  transition  from  forced  to  free  convection 
occurs  in  the  range  -0.02  <  Ri  <.  -0.05  and  more  particularly  at  Rl  =  -0.0317; 
and  (3)  that  the  interpolation  formulae  for  the  diabetic  influence  function  span 
the  transition  from  forced  to  free  convection.  Since  none  of  the  interpolation 
schemes  are  universally  accepted  and  verification  of  the  solution  for  K({/kn  by 
direct  observation  is  extremely  difficult,  the  empirical  results  of  this  study 
perhaps  do  not  represent  the  actual  ratio  of  the  exchange  coefficients. 
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